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We calculate the effective charge for multimagnon infrared (IR) absorption assisted by phonons in the parent insulating 
compounds of cuprate superconductors and the spectra for two-magnon absorption using interacting spin-wave theory. Recent 
measured bands in the mid IR [Perkins et al. Phys. Rev. Lett. 71 1621 (1993)] are interpreted as involving one phonon plus 
a two-magnon virtual bound state, and one phonon plus higher multimagnon absorption processes. The virtual bound state 
consists of a narrow resonance occurring when the magnon pair has total momentum close to (n,Q). 
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The discovery of high-Tc superconductivity in doped 
cuprate materials 0] has triggered a lot of attention in 
the parent quasi- two-dimensional spin-^ quantum anti- 
ferromagnets. So far complementary information on the 
antiferromagnetism has come from different probes like 
neutron scattering and Raman light scattering In 
this work we reinterpret recent infrared absorption mea- 
surements [|| (Fig. 0) in terms of phonon assisted multi- 
magnon absorption. The narrow primary peak in Fig. |l| 
is explain in terms of a long lived virtual bound state of 
two magnons here referred to as a bimagnon. These new 
states are narrow resonances occuring when the magnon 
pair has total momentum close to (7r, 0) and have a rezon- 
ably well defined energy and momentum in a substantial 
portion of the Brillouin zone. 

In principle IR absorption of magnons is not allowed 
in the tetragonal structure of cuprate materials. This is 
because in a typical two-magnon excitation the presence 
of a center of inversion inhibits any asymmetric displace- 
ment of charge and hence the associated dipole moment 
vanishes. However the situation changes if phonons are 
taken into account. In a process in which one phonon and 
two magnons are absorbed, the symmetry of the lattice 
is effectively lower, and the process is allowed. A simi- 
lar theory was put forward by Mizuno and Koide j| to 
explain magnetically related IR absorption bands found 
many years ago by Newman and Chrenko in NiO H . To 
the best of our knowledge, we present the first explicit 
calculation of coupling constant for phonon-assisted ab- 
sorption of light by multimagnon excitations and the line 
shape for two-magnon absorption. We consider a Cu-0 
layer but generalization of our results to other magneti- 
cally ordered insulators is trivial. 

Consider a three-band-Peierls-Hubbard model Q in 
the presence of an electric field (E) and in which for 
simplicity Cu atoms are kept fixed and O ions are al- 
lowed to move with displacements u i+s / 2 . Here i labels 



Cu sites and 5 = x,y, so that i + 8/2 labels O sites. 

Holes have an on-site interaction U d on Cu and U p on 
O, a Cu-0 repulsion U p d, on-site energies on Cu Ed, and 
on O £ p , and Cu-0 hopping t (0-0 hopping is neglected 
here). We define A = E p - E d + U pd , e = 2(E P ~ E d ) + 
U p . When an O ion moves in the direction of a Cu with 
displacement |u|, the corresponding on-site energy of Cu 
changes to first order by — (5\u\ and the corresponding 
Cu-0 hopping by a\u\. Opposite signs apply when the 
O moves in the opposite direction j^]. 

To calculate the coupling constants of light with one- 
phonon-multimagnon processes we first obtain a low en- 
ergy Hamiltonian as a perturbation expansion valid when 
t << A, e, Ud and when the phonon field and the electric 
field vary slowly with respect to typical gap frequencies, 



H = 



'J{E,{u i+s/2 })B, 
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Here -Bi+5/2 = SiSi+s with Si the spin operators, Hph 
is the phonon Hamiltonian containing spring constants 
and masses for the O ions (M) and Pph is the phonon 
dipole moment. The first term in Eq. p) contains spin- 
dependent fourth order corretions in t whereas fourth, 
second and zero order spin-independent processes are col- 
lected in the last two terms. As usual we calculate the 
superexchange, J, as the energy difference between the 
singlet and triplet states of the spins located at Cul and 
Cur in Fig. |^. We only need to consider the three con- 
figurations (A,B,C) of the L-R bond and E. Next we 
Taylor expand J to first order in E and {1*1+5/2} @i 

J = J + r;(u L - ii R ) - E[qiu + Aq A (2u - «l - %)] 

(2) 

Here A = 1 for configuration A and A = for configu- 
rations B and C. In each configuration the displacement 
of the central O and the electric field are parallel, i.e. 
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E = Ee, Uq = Moe. The direction of e is the same 
as the arrows at the bottom of Fig. ||. ul and Ur are 
only relevant in configuration A. ul = uli + «L2 — Uh3, 
Mr = — uri + ur2 + mr3- The numbering and the di- 
rection of the displacements are shown in Fig. [| The 
first term in Eq. (Q) is the superexchange in the absence 
of the electric and phonon fields, Jo = %?[tj^ + §]• The 
remaining quantities are a magnon-phonon coupling con- 
stant, r\ = = gp- /^[^(jj- + |) + ^] an d effective charges 
associated with one phonon and multimagnon processes, 
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a p ,i is the Cu-0 distance. Within a point charge estima- 
tion the parameter (3a p d ~ ^U p d- The dipole moment is 
obtained as P = — §j| and using Eq. (||) in the relevant 
configurations. We get up to fourth order in t, 



P = P 



lph 



lph+mag 



(5) 



The first term describes conventional phonon absorption. 
We define SB i+S / 2 = B i+S / 2 - (B i+S / 2 ) and its Fourier 
transform, SB p and the Fourier transform of u i+ g/ 2 , u p- 
After Fourier transforming, the dipole moment for one 
phonon and multimagnon processes for an in-plane field 
in the x direction has the form, 



Px, lph+mag = N[qiY^ SB s p u s XI 



p,6 



+A4g A ^sin(^)sm(^)^uf p ] 



(6) 
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and A = f . The case of an electric field perpendicular to 
the plane is obtained by putting A = and replacing u x 
by u s z . N is the number of unit cells. The first term is 
isotropic being present in any configuration. Looking at 
the cluster in Fig. |^ it can be understood as a spin depen- 
dent correction to the charge in Oo- Its physical origin is 
that fourth order corrections to the charges involve spin 
dependent processes. For example if the spins in Cul 
and Cur are parallel they cannot both transfer to Oo 
whereas if they are antiparallel they can. Fig. |(a) illus- 
trates a typical process efficient in configuration B. The 
second term is anisotropic being present for an in-plane 
field only. It originates from a "charged phonon" like 
effect fllO|] . Consider the configuration in which the elec- 
tric field and the displacement of Oo are both parallel to 
the Cul-Cur bond (A in Fig. ||), and a phonon in which 
the O's around Cul breathe in and the O's around Cur 
breathe out. (We don't need to consider zero momentum 
phonons to couple to light since it is the total momen- 
tum, magnons plus phonons which has to add to zero.) 
The Madelung potential in Cur decreases, and in Cul it 



increases, creating a displacement of charge from left to 
right that contributes to the dipole moment. Fig. ||(b) 
illustrates a typical process. But again this effect is spin 
dependent since if the two spins are parallel they cannot 
both transfer to Cur. 

The real part of the optical conductivity due to these 
processes is given by the dipole-moment-dipole-moment 
correlation function. Assuming, for simplicity, that only 
the wf/ p 's with the same S and 5' mix, and decoupling 
the phonon system from the magnetic system which is 
valid in lowest order in the magnon-phonon coupling we 
get (h = 1), 
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Here ui\\ p is the frequency of the u% p and phonons and 
Lu±p is the frequency of the u v xp and Uy p phonons. £j|| p 
can be associated with the frequency of Cu-0 stretch- 
ing mode phonons and lu± p with that of Cu-0 bending 
mode phonons. The supraindex in the Green functions 
indicates that the poles should be shifted by that amount. 
The absorption coefficient is obtain assuming weak ab- 
sorption as a = e ^L- <7 with e% the real part of the dielec- 
tric constant |J(b). 

To compute the magnon-magnon Green functions we 
use interacting spin-wave theory jl]J with a Holstein- 
Primakoff transformation. On the A (B) sublattice we 

put Sf(S^) = \J 25(1 — b\bi)bi and the corresponding 
Hermitian conjugates. 5 = 1/2 in our case and bi is 
a boson operator. With this definition the Hamilto- 
nian is invariant under the exchange of the sublattices 
and we don't need to distinguish between them. Ac- 
cordingly we work in the non-magnetic Brillouin zone. 
The Heisenberg Hamiltonian in this representation is 
expanded to zero order in 1/5 and normal ordered 
with respect to the non-interacting spin-wave ground 
state. The non-interacting part is diagonalized by 
the Bogoliubov transformation Q\ — Ukb\ — Vkb-k, 
where bk is the Fourier transform of bi and Uk = 



-sig(7fc)\/(l - Wfc)/(2w fc ), with 



+ U>k)/(2(J k ), V k 

7fc = 2/zJ2s cos(kd), ujk = \/l — 7fc and z is the co- 
ordination number (z — 4 in our case). We define 

9 Pqiq2 = ((Q±p+ qi Q± P - qi ;Q\ p+q2 Q\ p - q2 ))- Next we 

apply a standard RPA like decoupling for the equation 
of motion of g which can be solved for integrals over the 
Brillouin zone of g weighted with products of u's and v's 
Pi . Finally ((6B^_ p ; SB p )) can be written as a sum of 
such quantities with proper weighting factors The 
dashed-doted line in Fig. |l| gives the contribution to the 
line shape from magnon pairs with p = {tt, 0) without 
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the approximations that follow. A very sharp resonance 
occurs there indicating that a virtual bound state (bi- 
magnon) is formed. 

Analytic expressions for the Green functions can be 
obtained neglecting the contributions involving u's. The 
resulting error is small because at short wave length v is 
small to start with and at long wave lenght (low energies) 
the spectral weight is small |0|. This approximation 
does not shift appreciably the peak at p = (tt,0) and only 
decreases slightly the intensity Bj. The Green function 
takes the form ((5B x _ p ,6B*)) = j^G xx with, 



(~x TT — 



r (0) , 9 t r r (0) r (0) (r(°)\2] 
*-Jxx t ^oq^jxx '■Jyy ~ v^xy ) J 



1 + 'ZJoiGxx + Gyy ) + 4:Jq[Gxx Gyy — {Gxy ) 2 ] 



Where, 
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u l P + q u \ p - q cos (<l6)cos(qs') 
u - Ei p+q - Ei p _ q 



(8) 
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Ek = E m uj k , E m — zSJ (l + C/2S) and ( is the Oguchi 
correction ( w 0.158 in our case. Eq. (||) takes the 

In this 
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familiar form at p — (tt.O), G tt — 

case the imaginary part of G x °l starts to be different 
from zero at u> = B m (~ 2 Jo) and rises very slowly as 
the frequency is increased so that it is very small at the 
bimagnon energy resulting in a long life time or narrow 
peak. In the inset of Fig. |l| we show the imaginary part 
of the Green function from Eq. (|J) for different values 
of p. This shows how the bimagnon disperses around 
p = (tt, 0) it goes upwards on going towards (0, 0) and 
downwards on going towards (tt, 7r). This indicates that 
(-7T, 0) is a saddle point and hence it should give a Van 
Hove singularity when integrated over p [Eq. M)]. The 
dashed curve in Fig. [l] is the theoretical line shape in the 
same approximation and assuming that the anisotropic 
processes dominate (qi — 0) and treating the phonons 
as Einstein like with un — 0.08eV 13 1. It gives a sur- 
prisingly good fit for the primary peak. Such a good fit, 
especially for the width was not possible within RPA in 
the Raman case [pj . This can be partially reconciled by 
the fact that a structure that is artificially broadened at 
p = (7r, 0) (but still much narrower than the integrated 
line shape) does not change significantly the final result. 

Because of the Van Hove singularity the position of the 
p = (tt, 0) bimagnon peak coincides with the peak in the 
line shape and is given by 1.179i? m + w\\ = 2.731Jo + oj\\ ■ 
This provides an alternative way to estimate Jo- Wc 
found Jo = 0.121eV which is in good agreement with 
other estimates [p|JT^]. 

The phonon is very anomalous in orthorhombic 



La2Cu04 since it splits due to anharmonicitics 13 its 



partner being at u>t =0.06eV at room temperature. Pre- 
sumably this produces the shoulder observed at lower en- 
ergies in the experiments (Fig. |l]), although the distance 
to the primary is larger than expected which may be due 
to other phonons involved. This feature was assigned to 
direct two-magnon absorption [0 made weakly allowed 
by the lower lattice symmetry according to the results of 
Ref. [p| . However we found that the dipole moment for 
this process is directed in the direction bisecting an an- 
gle made by the Cu-O-Cu bond jl6| and hence can only 
contribute for a field perpendicular to the plane. 

The oscillator strength in Fig. |l] was adjusted to fit 
the experiments. A rough estimate of the strength using 
qa = -le and t\ = 5 gives a value ~ 4 times smaller than 
observed which is quite reasonable given the uncertainties 
involved. 

We interpret the side bands at higher energy to be due 
to higher multimagnon processes, neglected in the above 
approximations. In order to check both the validity of 
our approximations and the origin of the side bands we 
have also computed the absorption using exact diagonal- 
ization of a small cluster ||. The exact result confirms 
the Green function calculation for the two-magnon peak 
and shows also side bands corresponding to higher mul- 
timagnon processes which we associate with the higher 
energy side bands observed in the experiments. The rel- 
ative weight of the side bands seem to be smaller than 
in the experiments presumably because of finite site ef- 
fects or the presence of other processes in the magnetic 
Hamiltonian as was suggested in the Raman case [[lTj . 

For an electric field polarized perpendicular to the 
plane only phonons perpendicular to the Cu-0 bond con- 
tribute. We estimate the absorption to be roughly a fac- 
tor of 8 smaller than the in-plane contribution. From the 
experimental side the anisotropy seems to be larger. One 
should be aware that the cuprates are in a regime where 
covalency is not small with respect to typical gap ener- 
gies and hence a perturbation in t is helpful to identify 
the important processes and discuss trends, but quanti- 
tative estimations are to be taken with care . As higher 
orders in t are included we expect that the anisotropic 
contributions grow with respect to the isotropic ones. For 
example the charged phonon effects of Fig. 3(b) can be- 
come very efficient if the second hole forms a Zhang-Rice 
singlet with the hole already present in Cul since that 
process involves a much smaller gap. Note also that the 
larger the order in t the longer the range of the processes 
that contribute to the anisotropic charges whereas only 
local processes contribute to the isotropic charge. Longer 
range processes have in general very large form factors. 
These effects should give a stronger anisotropy in accor- 
dance with the experiments. This also partially justifies 
our simplifying assumption of taking q\ = 0. 

To conclude we have computed effective coupling con- 
stants of light with multimagnon excitations assisted by 
phonons and determined the line shape of the primary 
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peak. Our results explain recent measured absorption 
bands in the mid IR of parent cuprate superconduc- 
tors and demonstrate the existence of very sharp virtual 
bound states of magnons. 
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FIG. 1. Experimental data from Ref. || (solid line) and 
theoretical line shape for two-magnon absorption (dashed 
line) in La2Cu04. The dashed doted line is the contribu- 
tion to the line shape from the bimagnon at p — (it,0). The 
inset shows a blow-up of p — -Im(C„) at the energy of the 
bimagnon for different values of the total momentum. 

FIG. 2. Schematic representation of the cluster used in the 
calculations. Full dots represent Cu's and open dots O's. 
Thick arrows represent the spin, thin short arrows represent 
lattice displacements and thin long arrows represent the di- 
rection of the electric field. We have represented uo in A con- 
figuration. In general its direction is equal to the direction of 
the electric field. 

FIG. 3. Typical processes contributing to the isotropic (a) 
and anisotropic (b) effective charges. The meaning of the 
symbols is the same as in Fig. 0. 
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